AN ALGORITHM FOR FINDING LOW DEGREE RATIONAL 
SOLUTIONS TO THE SCHUR COEFFICIENT PROBLEM 
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Abstract. We present an algorithm producing all rational functions / with 
prescribed n + 1 Taylor coefficients at the origin and such that 1 1 / 1 1 oo < 1 and 
deg / < k for every fixed k > n. The case where k < n is also discussed. 



1. Introduction 

Let H °° be the Banach space of bounded analytic functions on the open unit 
disk ID) with norm ||/||oo := sup 2gD |/(z)| < oo. The closed unit ball S of H°° 
(sometimes called the Schur class) thus consists of analytic functions mapping D 
into its closure. The classical Schur problem which we will denote by SP„ consists 
of finding / G S having prescribed n + 1 Taylor coefficients at the origin. 

SP„: Given cq, . . . ,c n G C, find all functions f € S of the form 

f(z) = c + cyz + . . . + c n z n + 0(z n+1 ). (1.1) 

The problem has a solution if and only if the Pick matrix of the problem given by 



Pn 



T(c n , . . . , co)T(c r , 



,co) 



is positive semidefinite. Here and in what follows, / denotes the identity matrix 
of the size always clear from the context, and T(cq, . . . ,c n ) stands for the lower 
triangular Toeplitz matrix with the bottom row entries indicated in the parentheses: 



T(c„, . . . , c ) 



CO 
Cl 





CO 








CO 



(1.2) 



If Pn ^ is singular, then the problem SP n has a unique solution which is a finite 
Blaschke product of degree equal to the rank of P n . In what follows, we assume 
that the data set {co, . . . , c„} is such that P n > and we will call such a data 
set admissible. For an admissible data set, the parametrization of the solution set 
of the problem SP„ was established in [7] via the famous Schur algorithm which 
we now recall. Starting with Co, • • • , c n , define the numbers cj^ (J = 1, . . . , n; k = 
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0, . . . , n — j) from the following recursion: 



[ 4 0) " 




CO 
Cl 


and 


r -W+i) ■ 

ii+i) 
c i 


= 


c l 

Ji) 

c 2 


c (0) 




c n 








r (i) 



(i > 0), (1.3) 



where the matrix 



7 0")„0") 



3 0') M 



,1 



=0°° I 9 



is defined via formula (11.21). Let 



7j 



for i = 0,...,n. (1.4) 

If Co, . . . , Cn are the Taylor coefficients of an / G S, then the numbers ji constructed 
above are the n + 1 first Schur parameters of / and condition P n > is equivalent 
to |7j| < 1 for i = 0, . . . , n. The Schur algorithm relies on the following fact: 

A function f belongs to S and satisfies (jl.ip if and only if it is of the form 

z h( z ) + Co 



/(*) 



zcofi(z) + 1 



(1.5) 



for some f\ 6 S such that fi(z) = + z + 



+ c^I-lZ*- 1 + 0(z n ) where 



1 are the numbers defined via ([1.3 



Starting with a function /o := / € <S of the form {TUj) and applying recursion 
(jl.5p n times one gets a sequence of Schur class functions satisfying 



zfj+i(z)+c l 



(J) 



z f] + li z )+l0 



zc ( j) f j+ i(z) + 1 zjjf j+ i(z) + l 



U = o, 



and such that fj(z) = c ( 



0) 



-.0) 



0(« 



.,n) (1.6) 
) where c^' are 



the numbers defined via (| 1 . 3[) . Upon taking the superposition of linear fractional 
transformations (|1.6|) one gets the linear fractional formula 

which parametrizes all solutions to the SP„ where the free parameter £ := /„ runs 



through iS and the coefficient matrix & = 



A 
C 



D 
D 



is given by 



Q(z) = W Q (z)Wi{z)---W n {z) where Wj(z) 



z 



7j 
1 



(1.8) 



Motivated by engineering applications (where it is desirable for the solution / of an 
interpolation problem to be rational and of small McMillan degree), the rational 
coefficient interpolation problem (as well as its multi-point analogs) was considered 
in pQ with an additional constraint on the degree (complexity) of rational inter- 
polants. In what follows, the polynomials Nf and Df will denote the numerator 
and the denominator from the coprime representation / = Nf/Df of a rational 
function /. By deg/ = maxjdeg A^/,degD/} we mean the McMillan degree of /. 
The algebra of rational functions will be denoted by 1Z and we will let 

K k := {/ G K : deg / = k} and K< k :={f EK: deg / < k} . 
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Being adapted to the single-point case, the problem formulated in pQ is: 
RP n! fc: Given cq, . . . , c„ 6 C and k > 0, /md a£/ / £ lZ< k of the form 
The problem was solved in [T] and in [5] (for the matrix-valued case) as follows. 

Theorem 1.1. Let q denote the rank of the Hankel matrix H = [ci+j-i]; j >:L 
constructed from the given numbers Cj ( the matrix H is x if n is odd or 
x t| if n is even). Then 

(1) There is no f € 7?.fc satisfying /or every k < q or q < k < n — q. 

(2) There exists at most one function f of complexity k — q subject to (|1.1[) . 

(3) For every k > n — q, there are infinitely many solutions of the problem 
RP„.fc which are parametrized by the formula 

/ = T a .g:=- — — , 1.9) 

2l2lff + 2122 

where the coefficients 2ty are polynomials explicitly constructed from the 
data set and such that 

Stia' = 

a 22 



deg 



B21 



q and deg 



1 



and where the parameter g = N g /D g £ 1Z is such that 
degN g <k-q, degD g <k + q-n-l, 2t 2 i(Q)iV g (0) + 2t 22 (0)I? a (0) ^ 0. 

We refer to [2] for more details. In what follows, we use notation 

sn = snn, sn k = snn k and sn< k = sn n< k 

for the classes of functions in 1Z, lZ k and lZ< k respectively, which are bounded by one 
in modulus on ID). Upon imposing both i?°°-norm and complexity constraints (i.e., 
upon combining problems SP n and RP ni fe) we arrive at the following interpolation 
problem. 

RSP ni fe: Given an admissible data set co,...,c„ and k > 0, find all functions 
f S SK< k of the form (JTTJ). 

One may try to treat the latter problem using either formula (|1.9p or \1.7\ . In 
the first case, the complexity of / is completely controlled by the complexity of the 
corresponding parameter g and it suffices to pick up all parameters g with deg g < 
k — q leading via formula (| 1 .9|) to Schur-class functions /. However, this task is hard, 
since formula (|1.9|) does not control ||Tgt[<7] ||oo in terms of ||(/||oo- It may happen 
that a Schur class parameter g produces f ^ S and on the other hand, a Schur 
class function / G SlZ< k may arise from a non-Schur class parameter g. Although 
Theorem 11.11 guarantees that there are infinitely many functions / G H n +i-q of 
the form (|l.ip . it is not known whether or not one of them is of the Schur class. 
The question about the minimal possible k for which the problem RSP nj /j has a 
solution, is still open. 

It is not even clear from (|1.9[) that the problem RSP„.fc has solutions for k large 
enough. On the other hand, the affirmative answer for the latter question is read- 
ily seen from parametrization formula (|1.7p which in contrast to (|1.9j) . perfectly 
controls the 7J°°-norm of /: all Schur-class rational solutions to the problem SP„ 
arise via formula (|1.7p from some Schur-class rational parameter £. The complexi- 
ties of interpolants are controlled here to some extent. A straightforward induction 
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argument deduces from (|1.8|) that the coefficients A, B, C and D in ljl.7|) are poly- 
nomials of respective degrees deg A = n+ 1, dcgB < n, degC < n+ 1, dcgZ) < n 
and therefore, 

degT e [£] <n + l + deg£. (1.10) 

Letting £ in (|1.7| to run through the class of constant functions (not exceeding one 
in modulus), one gets a family of solutions / of the problem RSPn,n+i, but not all 
the solutions. It turns out that zero cancellations may occur in (|1.7p due to which 
some solutions to the RSP„ Jl+ i may arise from non-constant parameters. We also 
observe that the parameter £ = leads via (|1.7p to the function Te[0] = B/D £ 
SlZ< n which is therefore, a solution to the problem RSP n>n . The next example 
shows that this function might be the only solution to the RSP n! „. 

Example 1.2. Let |co| < 1 and Cj = for j = 1, . . . , n. With this data, the 
problem RSP„ „ has only one solution / = cq. This follows from Theorem 11.11 
since in this case q = 0. 



Otherwise (that is, if Cj at least for one j > 1 so that q > 1), Theorem ll.il 
guarantees the existence of infinitely many functions / £ TZ< n of the form (jl.ip . 
at least one of which (Te[0]) belongs to SlZ< n . As was shown in [4]- [6], the set of 
such functions is infinite and can be parametrized by polynomials a with deg a < n 
and with all the roots outside B. More precisely, for every such a, there exists a 
unique (up to a common unimodular constant factor) pair of polynomials a(z) and 
b(z), each of degree at most n and such that 

(1) |a(z)| 2 - |6(z)| 2 = |a(z)| 2 for \z\ = 1 and 

(2) the function / = b/a (which belongs to SlZ n by part (1)) satisfies (jl.ip and 
therefore, solves the RSP„ in . 

The objective of this note is to present an alternative parametrization of the solu- 
tion set of the problem RSP nj fc (see Theorem 11.31 below) which relies entirely on 
parametrization formula (|1.7|) . Some elementary analysis of the Schur algorithm 
will relate complexities of £ and degTe[£] more accurately than in (| 1 . 10[) : this in 
turn, will allow us to describe all parameters £ G S1Z leading via formula ljl.7|> to 
solutions / of the problem RSP Wi fc (these parameters will be called admissible). 
Explicit construction of these parameters is given below in terms of certain algo- 
rithm which seems to be quite efficient and simple from the computational point of 
view. Here is the Algorithm: 



Step 1: Given co,...,c n , compute the numbers 70, 7i, • • • , 7n by formula ([1.4 
using iteration (|1.3p . 



Step 2: Using the numbers 70, . . . ,7„ compute the polynomials 

n n 

A n (z) =^2a j z 3 and B n (z) = ^& J z : ' (1-11) 

j=0 j=0 

from the system of recursions 

A (z)=j n , B Q (z) = l, 
A j+1 {z)=zA j {z)+*y n - j - 1 B j (z), . , (1.12) 

It is readily seen that Bj(0) = 1 for j = 0, . . . , n. In particular, bo = B n (0) = 1. 
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Step 3: Using the coefficients a,j, bj from define the lower triangular Toeplitz 

matrices 



A = 



a n 







1 

h 



(1.13) 



a\ a 2 ■ ■ ■ a n 

and compute the lower triangular Toeplitz matrix 

R = T{r u r 2 ...,r n ) := B X A. (1.14) 

The three first steps are preliminary and can be carried out in finitely many steps. 
The last step tells which parameters £ in (|1.7[) should be taken to get solutions to 
the problem RSP,,^.. We first consider the case where k = n. 

Step 4: For any n-tuple {ax, ...,a„} of complex numbers, compute the function 

00 + [3 lZ + . . . + fin^z"- 1 



£{z) 



-»0 

CtQ + Oi\Z 



a n z' 



where ■ ■ • , n -i are defined by 



Pn-1 








= R 




00 







(1.15) 



(1.16) 



where R is given in (|1.14[) and ao is such that 

The main result of the paper is the following theorem; the proof will be given in 
Section 2. 



Theorem 1.3. Let £ be constructed as in Step 4 and let be as in (|1.8[) . Then the 
function f = Tq[£] (|1.7p solves the problem RSP„ jn and conversely, all solutions 
of the RSP„.„ arise in this way. 

Remark 1.4. The only relatively uncertain part in Step 4 is the choice of ao- 

n 

However, it is readily seen that for any ao satisfying |ao| > ^~^(|aj| + the 

i=l 

function £ in (|1.16p belongs to the Schur class which immediately gives infinitely 
many solutions of the problem RSP„ in . To be more precise, let us write (|1.15l) as 

P{z) 



£(z) 



ao + zQ( Z y 



Oil 



and let 



where P(z) = O + 1 z + . . . + p n ^ 1 z n ~ 1 and Q(z) 
D(c, r) denote the disk of radius r centered at c. Then the set of all admissible 
ao's (for already chosen a±, . . . , a„ and /3o, ■ ■ ■ , Pn-i) is the exterior (complement) 
of the domain 17 defined as 



n 



in 

\z\<l 



K-zQ(z),\P(z)\). 



Remark 1.5. It follows from (|1.15[) that a parameter £ leading to a solution of 
the RSP„ jn has to satisfy £(oo) = 0. Thus, £ = is the only admissible constant 
parameter for the problem RSP n! „ . Combining this fact with (|1.10[) . we conclude 
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that every other constant function £ £ S leads via (11. 7|) to a solution of RSP 



n,n+l ■ 



As we have already seen, in contrast to the case n = k, the existence of infin- 
itely many solutions of the problem RSP^fc with k > n is immediate. However, 
the description of all solutions is even somewhat more complicated. We get this 
description by an appropriate modification of Step 4 as follows. 

Step 4': Let k > n be fixed and let and R be as above. All solutions f to the 
problem RSP,^ are obtained via formula (jl.7p where the parameter £ is either any 
function from S1Z<k-n-i or a function from SlZ< k of the form 

c , x Pn-k + Pn-k+lZ + . . . + f3 n -lZ k ~ 1 

£(z) = — (1.17) 

a n -k + a-n-k+iz + . . . + a n z K 

where the coefficients a n -h+i, ct n -k+2, ■ ■ ■ ,ot n and f3 n -k, Pn-k+i, ■ ■ ■ , P-i are picked 
up arbitrarily, after which the coefficients (3q, ■ ■ . , (3 n -i are defined as in (|1.16[) and 
where after all, the coefficient a n -k is chosen so that the function £ of the form 
(|1.17|) belongs to the Schur class S. 

Justification of Step 4' will be given in Section 2. In Section 3 we will present a 
version of Step 4 suitably modified for the case where k < n. There we will explain 
the reasons (by means of parametrization formula (|1.7p ) for which the algorithm is 
not efficient for k < n. 



2. Proof of Theorem 11.31 

In this section we justify the algorithm presented in the previous section. Let 

e fc (*) := W n - k {z)---W n {z) (2.1) 

where the factors Wj are defined in (|1.8|) . Comparing (|2.1[) and (|1.8p we see that 
n equals the coefficient matrix of the transformation (|1.7p . It is not hard to 
check by induction that 0^ is of the form 



Qk{z) = 



zBl(z) A k (z) 
zA{{z) B k {z) 



(2.2) 



where the polynomials A k and B k are constructed from system (|1.12|) and where 
Al and B\. are defined as follows: 



A\{z) = z k A k (l/z), B{{z) = z k B k (l/z). (2.3) 
• — A -' 

_ D e 

0„ = into fL7) 



Let us take any £ = jf- 6 STZ and substitute it together with formula (|2.2[) for 



zBl(z)N £ (z) + A n (z)D £ (z) 
' [! zAi{z)Ns{z) + B n {z)D s {z)' 



Remark 2.1. The numerator and the denominator in (|2.4[) do not have common 
zeros and thus, 

N f = zBl N e + A n D e and D f = zA^Ng + B n D s . (2.5) 
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Proof: Taking determinants in (11. 8|) . (12. 1|) and ()2.2jl (with k = n) gives 

S„(z)^(z)-A„(z)4(z) = i-det6„(z) (2.6) 

1 n n 

= -■l[detW j (z)=z n .l[(l-\ 7j \ 2 ). 

Z j=o j=o 

Therefore, the only possible common zero for the numerator and the denominator 
in (f2Tf is z = 0. But if this is the case, we then have B n (0)D e (0) = D £ (0) = 
which is impossible since the Schur function £ cannot have a pole at the origin. □ 

We shall now compare McMillan degrees of / and f\ in formula (| 1 . 5[) . 

Lemma 2.2. Let f e SK be of the form ([Q]) . Then deg / - 1 < deg A < deg/. 

Moreover, 

deg/i=deg/ /i(oo) = /(oo)^l/co (2.7) 

and 

deg/ 1 = deg/-l ^=> /(oo) = !/£*,. (2.8) 



Proof: Take /i in the form f\ — Nf 1 /Df 1 and rewrite (|1.5|) as 

z^^+cp^^z) _ F(z) 
7W +^0*) G(z) V»> 

from which we see that deg Nf < deg fi + 1, deg Df < deg /i + 1 and thus, deg / < 
deg /i + 1 . Now let us take / in the form / = -^p and solve equation (|1.5| for /i : 

= (iy /( z) -coD^z))/, (210) 

m; D f (z)-c N f (z) { ' 

Since Co = /(0) = Nf(0)/Df(0), it follows that the numerator in (|2.10p is a polyno- 
mial of degree not exceeding deg / — 1 . Therefore, deg-ZV^ < deg/, deg Df ± < deg / 
and thus, deg/i < deg/. This completes the proof of the first statement. 

Since there are only two possibilities for the value of (deg / — deg f\ ) , statements 
(|2.7[) are equivalent to (|2.8p . We next observe that the polynomials F and G in (|2.9|) 
do not have common zeros (the proof is the same as in Lemma 12.21) and therefore 
we can conclude from (12.91) that 



deg / = max{degF, degG}. (2.11) 

Now we verify (|2.7[) (or (|2.8p ) separately for the following three cases. 

Case 1: Let degZ?^ > degiV^ + 1. Then it follows from (j2l?|) that deg/ = 
degDfi = deg/i and on the other hand, /i(oo) = and /(oo) = Co ^ l/co- 

Case 2: Let degD h < degiV^ + 1. Then deg / = degN fl + 1 = deg/i + 1 and 
on the other hand, /i(oo) ^ and /(oo) = 1/cq. 

Case 3: Let deg Df 1 = deg iV^ + 1. Let ao and bo be the leading coefficients of the 
polynomials Nf 1 and D t x respectively. Then the leading coefficients of F and G are 
a + c &o and co^o + foo = 0, respectively. Assuming that degF < degNf 1 + 1 and 
deg G < deg Nf ± + 1 we have a + co^o = and coao + b = which gives a = b = 
which is a contradiction. Therefore, max{deg F, deg G} — deg iV^ +1 and by (|2.11jl , 
deg / = degTV/j + 1 = degD/j = deg/i. Finally, since deg-D/j = deg Nf 1 + 1, we 
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have /i(oo) = and it follows from (|2.9|) that f(oo) = — which is not equal 

to 1/co, since bo ^ and |cq| ^ 1. □ 



Let us apply the backward Schur algorithm (|1.6p to a function £ G iS7\Lfc by 
letting 

f n+1 =£ and /,■(*) = for j = 0,...,n. (2.12) 

Lemma 2.3. If deg fi = deg/i+i + 1, i/ien deg fj = deg/j+i + 1 /or every j < i. 
If fi(oo) = 0, i/ien fj(oo) = and degfj = deg fi /or every j > i. 

Proof: If deg/j = deg/ i+1 + 1, then by virtue of (|2.8p (with /, /i and c 
replaced respectively by fi, f i+ i and 7j) we have /i(oo) = i ^ 0. Then again by 



(applied to the new triple fi-i, ft and 7,_i) we get deg/i_i = deg fi + 1 and 
therefore, /i-i(oo) = 7^ 0. The first statement then follows by induction. 

We now assume that /,(oo) = 0. Since fi(oo) ^ =-, we conclude from (|2.7p that 
/i + i(oo) = and deg/i + i = deg fi. The induction argument completes the proof 
of the second statement. □ 

Proof of Theorem 11.31 Let / be a solution to the problem RSP„ in , i.e., / is a 
rational Schur-class function of degree at most n satisfying equality (jl.ip . Then / 
is of the form Q1.7|) for some rational Schur-class function £ or equivalently, f = fo 
is obtained from £ — f n via recursion (|2.12[) . Then we necessarily have 

deg£<n and £(00) = 0. (2.13) 



Indeed, deg / > deg £ by Lemma 12.21 and since deg / < n by the assumption, the 
first relation in (|2. 13|) follows. If we assume that £(00) 7^ 0, then we get by virtue of 
(|2.8|) , that deg f n -i — deg £ + 1 and then we also have deg / = deg £ +n + 1 > n + 1 
(by the first statement in Lemma I2.3[) which contradicts the assumption. Thus, 
£ (00) = 0. Due to (|2.13p we can take £ in the form (|1.15p . i.e., we can let 

n— 1 n 
j=0 j=0 

It remains to show that the coefficients <Xi and Pi are related as in (| 1 . 16|) . Observe, 
that the polynomials A n and B n constructed in (11.12p are of degree at most n; 
we take them in the form p. lip so that the reflected polynomials A\ and B\ (see 
(|2.3p ) are given by 

n n 

At{z) = J2<in-jZ j and B»(z) = £ b n ^z j . (2.15) 
3=0 j=0 

Substituting (fTTTTjl . (|2~T4| and (f2~T5|) into ((275]) we get 

n-l /n-l-1 \ 

fcO \ ]=0 ) 
n-l (n-l-\ \ 
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where Pi and P2 are polynomials of degree at most n. The two latter formulas 
imply that deg / < n if and only if 

n-l-l 

^ (b n -l-j-\f3n-j-i + ai+j+xoin-j) = (£ = 0, . . . ,n - 1), 

Jl— 1 

2J (S n -<e-i-i/3n-j-i + &<+i+i"n-i) = ° (^ = 0, ... 
3=0 

Making use of the Toeplitz matrices 

A = T(ai,a 2 , . . . ,a n ), B = T(bi,b 2 , . . .,&«), 

A = T(a„_i, a„_ 2 , . . . ,ao), B = T(6 rl _i, fe„_ 2 , ■ ■ ■ 
and of the vectors 

" Pn-l ~ 

and (3 = 



,n-l). 



(2.16) 
(2.17) 

(2.18) 



Oil 



0o 



(2.19) 



(2.20) 



one can write equations ()2.16j) and (|2.17|) in the matrix form as 

B/3 + Aa = and A/3 + Ba = 

respectively. Since 60 = B n (0) = 1, the matrix B is invertible. Then we get from 
the first equation in (|2.20p 

/3 = -B^Aa = -Ra (2.21) 

which is the same as (|1.16j) . We thus showed that every solution / to the problem 
RSP ?l! „ can be obtained via the Schur algorithm from a parameter £ G S of the 
form (fTT5l) . (flT6)) . 

To show that any such parameter is admissible, we have to verify that the vectors 
a and (3 related as in (I2.21|) satisfy both equations in (|2.20|) . The first equation 
is clearly equivalent to (|2.21|) . Substituting (|2.21j) into the second equation and 
taking into account that all the matrices in (|2.18p commute, we get 

-AB 1 Aa + Ba = r 1 (BB-AA) a. 



A/3 + Be* 



(2.22) 



We next substitute formulas (jl.lip and (|2.15[) into (|2.6p and examine the coefficients 



of 



for I = 0, . . . , n — 1 to get equalities 



{bn+j-ibj - a n+ j-£a,j) = (£ = 0, 



1), 



(2.23) 



which can be written in terms of matrices (|2.18p as BB = AA. We now conclude 
from [|2.22p that the second equation in (|2.20p is satisfied. Thus, for every £ S S 
of the form (fTTT5| . (fTTTB]l . the coefficients ecu fa satisfy equalities ([2~T6]) . (f2TT7|) 
(i.e., equalities (|2.20p ). which in turn guarantees that the McMillan degree of the 
function / obtained from £ via the Schur algorithm, does not exceed n. Since this 
/ belongs to S and satisfies (|1.1|) . it solves the problem RSP„ !ra . □ 

Justification of Step 4': Let k > n be a fixed integer. Every solution / to 
the problem RSP„,fe is of the form (|1.7p for some rational parameter £ G S1Z 
with deg£ < k. We have either £(00) ^ or £(00) = 0. In the first case, 
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deg / = deg £ + n + 1 (by Lemmas 12.21 and 12. 3j) and therefore, deg £ < k — n—1. On 
the other hand, for every £ E S7Z<k- n -i, it follows from (|1.10|) that degTe[£] < k. 
In the second case, we can take £ in the form (|1.17[) . that is to let 



N £ (z) 



E 

3=0 



Pn-k 



and Dz(z) 



k 

E 

3=0 



Substituting the latter formulas along with and (|2.15p into (|2.5p we get the 

formulas for Nf and Df as in the proof of Theorem 11.31 but with the factor z k+1 
(rather than z n+1 ) on the left and with polynomials P\ and P 2 of degree at most 
k. Then we conclude that deg/ < k if and only if conditions (|2 . 20[) hold which is 
equivalent to (|2.2ip . 

3. Concluding remarks 

In conclusion we present a version of the main algorithm for the case where 
k < n. The three first steps are the same as before; the last step describing all 
admissible parameters in parametrization formula (|1.7|) is the following. 

Step 4": Let k < n be fixed and let and R = T(ri , Ti , ■ ■ ■ , r n ) be as above. 
All solutions f to the problem RSP n! fc are obtained via formula (|1.7p where the 
parameter £ is a Schur- class function of the form 

N £ (z) 0o + p lZ + . . . + Pk-xz^ 1 



£{z) 



Ds(z) 



Q'O + OL\Z + . . . + ttfcZ* 



where the coefficients ctQ, . . . ,ctk and /3q, ■ ■ ■ , Pk-i satisfy the system 



' Pk-l ' 




Tn • 

r n -i r n 


■ 
• 




Cik 


Pi 






. 






Po 




1 n — k-\- 1 — k 


• r « 







ri 
r-2 



''2 

•'•3 



Tk+l 
Tk+2 



Vn—k Tn—k-\-l 



Oik 
(Xk-1 

a 



(3.1) 



(3.2) 



(3.3) 



(the matrix in p.3p is of Hankel structure). 

Proof: As in the proof of Theorem 11.31 we first observe that every solution 
/ of the problem RSP n ^ is of the form (jTTTJ) for some £ £ SlZ<k subject to 
£(00) = 0. Therefore, £ can be taken in the form (|3.1[) . Substituting (|1.11[) . 
(|2.15p and (|3.ip into (|2.5|) we now get Nf and Df the polynomials of degree at 
most n + k. Then equating the coefficients of z k+t of these polynomials to zero for 
j = £,..., n — 1, we get necessary and sufficient conditions (similar to (|2 . 1 6|) and 
(|2.17p ) for deg / = maxjdeg Nt, deg Df} not to exceed k. These conditions are 

-i,fe} 



min{n-£-l.fe-l} 

' E * 

3=0 

min{n-£-l,fc-l} 

E 

3=0 



'n-i-j-l 



Pk- 



3-1 



min{7i- 

E 



a e+ j +1 a>k-j — 0, 



3=0 

min{n — i— Ijfe} 



Q>n-i-j-iPk-j-l + 



E 

3=0 



^1+3 + lCtk-j 



= 
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(£ = 0, . . . , n) and it is not hard to see that they can be written in the matrix form 

flmni) as 

B/3 + Aa = and A/3 + Ba = (3.4) 

respectively where the matrices A, B, A and B are the same as in (|2.18p and where 
now 



oik 



ao 








and (3 



fc-i 

A) 
o 
o 

o 



(3.5) 



Since BB = AA, it follows as in the proof of Theorem II. 3[ that a. and (3 solve 
the system (|3.4[) if and only if they are related as in (|2.2ip . Substituting (|3.5|) 
into (|2.2ip and comparing the k top entries in the obtained equality, we get (|3 . 2[) ; 
comparison of the n — k bottom entries gives (|3.3| . □ 



Remark 3.1. Although Step 4" looks very similar to Step 4 in Section 1, in fact it 
is much less efficient. Let us demonstrate this by the case where k = n — 1. Then 
condition (13.31) takes the form 



nctn-i + r 2 a n -2 + ■ • • + r n _iai + r n a Q = 



(3.6) 



and if r n ^ 0, then ao is uniquely determined by ax, ... , a n _i. The problem is to 
describe all the tuples {or, . . . , a„_i} (which now are the only free parameters) for 
which the function 

/3 +/? 1 z + ...+/3„_ 2 z"- 2 



£{z) 



(3.7) 



a + ol\z + . . . + a n -iz n 1 

with the coefficients ao, Po, - ■ ■ determined by formulas (|3.6p and (|3.2p (with 

k = n — 1), belongs to the Schur class. The problem is hard; at the moment we 
even do not know necessary and sufficient conditions for the existence of at least 
one such tuple. 
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